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Abstract 

We present a novel method for studying the shape of the Z boson transverse mo- 
mentum distribution, Qt, at hadron colhders in pp/pp Z/j* — > l^l~. The Qt 
is decomposed into two orthogonal components; one transverse and the other par- 
allel to the di-lepton thrust axis. We show that the transverse component is al- 
most insensitive to the momentum resolution of the individual leptons and is thus 
more precisely determined on an event-by-event basis than the Qt- Furthermore, we 
demonstrate that a measurement of the distribution of this transverse component 
is substantially less sensitive to the dominant experimental systematics (resolution 
unfolding and Qt dependence of event selection efficiencies) reported in previous 
measurements of the Qt distribution. 
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1 Introduction 



The shape of the Z boson momentum distribution transverse to the beam 
direction {Qt) at a hadron collider tests the predictions of quantum chromo 
dynamics (QCD), since non-zero Qt is generated through radiation from the 
initial state partons. A good understanding of electroweak vector boson pro- 
duction is important in precision measurements (e.g. top quark and W boson 
mass) and in Higgs boson and new phenomena searches at hadron collider ex- 
periments. In many of these searches the signal and/or background processes 
involve electroweak vector bosons produced in association with jets. 
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At low Qt {Qt ^ Q), where the emission of multiple soft gluons is important, 
calculations in fixed order perturbative QCD diverge. There exist resummation 
techniques, in which singular contributions from all orders of are resummed 
to give a finite result. Resummation was first apphed to the Drell-Yan process 
by CoUins, Soper and Sterman (CSS) [1]. The resummation is carried out in 
impact parameter (6) space and includes a non-perturbative (NP) form factor 
that needs to be determined from data. Brock, Landry, Nadolsky and Yuan 
(BLNY) proposed the following parameterisation [2]: 
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where Xi and Xj arc the fractions of the hadron momenta carried by the initial 
state partons, Qq is an arbitrary scale set to 1.6 GeV and the parameters gi 
need to be fitted from data. 

Using the BLNY NP form factor, the CSS formalism was able to describe si- 
multaneously Tevatron Run 1 Z data and lower Drell-Yan data in a global 
fit of the parameters gi [2]. The Qt distribution at the Tevatron (Q^ ~ M|) 
is sensitive to g2 and insensitive to gi and g^. A larger value of g2 corresponds 
to a harder Qt distribution. The CSS formalism is implemented in the next 
to leading order (NLO) event generator ResBos [3]. In Run II the D0 Collab- 
oration reported a Qt measurement in the di-electron channel with 1 fb^^ of 
data [4]. For low Qt (Qt < 30 GeV), the D0 data is, within the measurement 
uncertainties, well described by the CSS/BLNY formalism. 

At low Qt the overall uncertainties were dominated by the parton distribution 
functions (PDFs) and the following experimental systematics [4]: 

• Unfolding the Qt measurement to account for the resolution in the mea- 
surement of the Et of the electrons. 

• Correcting for the Qt dependence of the overall event selection efficiency. 

As a result of these substantial experimental systematics, the low Qt region 
was not much better measured in this 1 fb~^ Run II analysis than in the 

100 pb~^ Run I analysis [5]. In both analyses, a measurement was made of 
g2. The measurements: 0.59 ± 0.06 GeV^ for Run I and 0.77 ± 0.06 GeV^ for 
Run II, have comparable uncertainties ^ . 

An observable that is sensitive to the Qt, but less sensitive to these experi- 

^ The D0 Run I and Run II measurements cannot be directly compared since they 
used different PDFs and the Run I measurement used the the Ladinsky-Yuan (LY) 
parameterisation [15] of the NP form factor as opposed to the BLNY parameteri- 
sation. Only the third term in gig^ is different and any shift in the fitted 52 is far 
smaller than the uncertainties on these measurements. 
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mental systematics would be beneficial. In this work, Monte Carlo studies of 
an optimal observable, a^, are presented. The observable has previously 
been used in the selection of I'l^vu final states at LEP by the OPAL collab- 
oration [6]. The UA2 Collaboration used a similar observable, p^, which we 
refer to as 6t, in a Qt measurement [7]. 



2 Constructing the Observable 

The measured Qt is highly sensitive to the lepton pt resolution. Our goal is to 
build an observable that is less sensitive to this resolution, whilst still sensitive 
to the Qt- We keep in mind the fact that collider detectors generally have far 
better angular resolution than calorimeter Et or track pt resolution. 




Fig. 1. A schematic representation in the transverse plane, of the construction of ay 
and ol in a typical leptonic Z decay. The hadronic recoil is expected to have equal 
and opposite transverse momentum to the Z. 

For events with di-lepton azimuthal separation, A0'^ > |, the Qt is decom- 
posed into orthogonal components as follows (See figure 1): 

• The thrust axis is defined as: t = io\ where pt'-*-' is the trans- 

verse momentum vector of lepton i. The two leptons have equal momentum 
transverse to this axis. 

• The transverse momentum vector of the di-lepton system, Qt = Pt^^^ +Pt^'^^ , 
is decomposed into components transverse to the axis, = \Qt x t|, and 
aligned with the axis, = Qt ■ t- 

For events with < |, is set equal to the Qt, while maintains the 
same definition for all values of A0'^ 

At low Qt, A0" ~ TT, hence the uncertainty on is approximately the 
uncertainty on the individual lepton pt^s multiplied by the sine of a small 
angle. In contrast, the uncertainty on (and thus also Qt) is approximately 
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the uncertainty on the individual lepton pxs multiphed by the cosine of a 
small angle. 

An alternative observable is discussed, br, whose construction is identical to 
that of qt except that the decomposition is performed relative to the di-lepton 

perpendicular bisector axis: b = ,^^~^'^J where r = |pt^^^|/|pt^^^|- The 

\PT —rpT I 

two leptons have equal acoplanarity with respect to this axis. In an event in 
which the leptons have equal px, the values of ax and bx are equal. No study 
is presented of the component longitudinal to the axis, 6^. 

As discussed below, the relative sensitivity to lepton px mis-measurement of 
ax and bx depends on the correlation between the lepton px and its resolution. 



3 Monte Carlo Simulation 



Monte Carlo (MC) events are generated using PYTHIA [8], which treats at 
leading order (LO) the process pp — > — > plus up to one jet, at a 

centre of mass energy of 1.96 TeV, and mass between 60 and 130 GeV. pythia 
uses additional parton showering to simulate "soft" transverse momentum 
generation. Although this study is carried out with Z IJ'~^IJ'~ events, the 
idea is applicable to studying the di-electron channel. 

In order to simulate the imperfect muon px resolution of a detector, Gaussian 
smearing of width 0.003 GeV"''' is applied in 1/px to both muons, which is 
approximately the design muon px resolution of the Run II D0 detector [9]. 
The design tracking resolution of other current and future hadron collider de- 
tectors; CDF, ATLAS and CMS [10-12] vary between 10"=^ and ^ 10"^ 
GeV~^ depending on the track px and pseudo-rapidity (rj). We refer to the 
constant S{l/px) form of the resolution px dependence as muon-like. An alter- 
native resolution dependence is studied; electron-like resolution with the form: 
5px/px = 0.15/p^/^, which would be the form expected for a calorimeter-based 
measurement appropriate for electrons. Gaussian smearing is applied to the 
lepton azimuthal angle of width 0.0005 rad., which is the typical resolution 
of a detector. 

In order to study the dependence of event selection efficiency on Qx, ar, ol 
and bxi the following simple cuts are applied to the event sample. These are 
representative of the cuts typically applied to select Z decays at a hadron 
collider. 

• The di-muon invariant mass must be between 70 and 110 GeV. 

• Kinematic cuts on both muons: px > 15 GeV and \ri\ < 2. 
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• An isolation cut on both muons: fiso < 2.5 GeV, where fiso = J2Pt{'^R < 0.2) 
and AR = ^ {AriY + (A^)^ is the radius of a cone around the candidate 
muon. The sum is over all particles excluding both the candidate muon and 
any neutrinos. 

PYTHIA is a LO event generator and is not expected to give a good description 
of the Qt distribution in real data. The events are re-weighted in Qt and Z 
rapidity (y) to match the prediction of ResBos using the default settings, which 
is in good agreement with D0 Run II data [4] in the region of low Qt {Qt 
< 40 GeV). Event samples are generated using ResBos, and the re- weighting 
procedure is carried out as follows: 

• Two dimensional histograms in Qt and \y \ are produced for both the pythia 
and ResBos event samples. 

• Dividing the ResBos histogram by the pythia histogram gives an event- 
weight histogram. 

• The PYTHIA events are now given an event-weight based on their Qt and 
\y\ such that the distributions of these variables are the same as the ResBos 
prediction. 

Hereafter, unless otherwise stated, MC refers to pythia re-weighted to Res- 
Bos. 



4 Sensitivity to Lepton pt Mis-measurement 
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Fig. 2. Mean resolution, \X^f.t — XgenM^gen for each of the observables, 
Qt-, clt, cll and 6t as a function of QTgen, with (a) muon-like resolution 
{6{1/pt) = 0.003 GeV"-*^) and (b) electron-like resolution {6pt/pt = 
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Two quantities are defined: the smeared or detector /eve/ quantity, Xdet, and 
the unsmeared or generator level quantity, Xgen, where X corresponds to Qt, 
gt, Ojl or hx- 

Figures 2 (a) and 2 (b) show the mean resolution, \Xii^t — Xg^^l/ ^gen as a 
function of Qxgen for muon-like and electron-like resolution respectively It can 
be seen that for low to moderate values of Qt {Qrgen < 50 GeV), qt and 6r 
are significantly better measured than or Qt- Either or Bt are therefore 
particularly well suited to studying Z production at low to moderate Qt- In 
the region of low Qt {QTgen < 15 GeV), gt and 6r have similar resolutions for 
either the (a) muon-like or (b) electron-like cases. At increasing Qt, ot is more 
suited for (a), and 6^ rnore suited for (b). The construction of ensures that 
the acoplanarity angle to the thrust axis is smaller for the larger pt lepton 
which tends to have poorer resolution in the muon-like case. In the following 
sections, for brevity, we discuss the case of muon-like smearing. 



5 Event Selection Efficiency Dependence 



Figure 3 shows the dependence of the event selection efficiency (separately 
for the cuts on muon \r]\, pt and isolation) on the generator level Qt, clt, gl 
and 6t- The cuts are apphed in the following order: \r]\; pt] isolation. The 
efficiency dependence for each cut is calculated having applied all previous 
cuts. For large Qt, the muons tend to be more central, increasing the rj cut 
efficiency. The same correlation is apparent for and ql, although weaker 
than for Qt- 

The muon pt cut dependence on is fiat in the range considered. Conversely 
large values of ai,, generate an asymmetry in the pxs of the two muons and 
tend to push the lower pt muon below the cut threshold. The dependence 
on the isolation cut is substantially flatter for than for ql- Large ax, corre- 
sponds to a high Pt hadronic recoil with a large fraction of its pt aligned along 
the thrust axis and thus possibly lying within the isolation cone of one of the 
two muons. There is no such dependence on ot, since or is the component of 
the recoil pt transverse to the thrust axis. 

In summary, the efficiencies of the cuts on muon \ri\, px and isolation depend 
less strongly on qt than Qt- The dependence of 6t on each of the cuts is 
similar to that of oy. 
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Fig. 3. The dependence of event selection efficiencies on the generator level (a) ot, 
(b) 6t, (c) aL and (d) Qt in Z ^ MC events. The efficiencies e^, Ep^ and 

eiso (evaluated relative to the previous cut in the order; |r/|, pT-, isolation) are shown 
separately for cuts on muon \'q\, pT and isolation respectively. Note that the range 
of Qt shown in this figure is ^/2x larger than for qt and a^. 

6 Systematic Uncertainties on the Unfolded and Efficiency Cor- 
rected Distributions 



Figure 4 compares (separately for Qt, aT, aL and 6t) the generator level 
distributions without selection cuts and the detector level distributions with 
selection cuts. The detector level Qt and ai distributions are substantially 
affected by the detector resolution and event selection efficiency. In contrast 
the distributions in aT and Bt are almost completely unaffected. 

In order to extract the underlying shape, the measured distributions would 
need to be unfolded for the detector resolution and corrected for the event 
selection efficiency with associated systematic uncertainties. These uncertain- 
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Fig. 4. Detector and generator level distributions for (a) ay, (b) 6t, (c) ql and (d) 
Qt- The detector level distributions are for Gaussian smearing in 1/pt of width 
0.003 GeV~^, and all selection cuts are applied. The generator level distributions 
do not include selection cuts. The lower halves of each plot show the fractional 
differences. 



ties would be substantially smaller for or hx than for Qj- or a^. In fact, for 
any realistic detector resolution, the measured or hr distribution describes 
the underlying distribution within a few percent without any unfolding or ef- 
ficiency correction at all. The statistical sensitivity of qt or 6^ to the shape 
of the underlying Qt distribution is also enhanced by the lack of resolution 
smearing compared to Qt- 
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Fig. 5. The generator level prediction of the ut distribution of pythia (re- weighted 
to ResBos) for two different g2 values. 

7 Uncertainties in Measurements of NP Phenomenological Param- 
eters 



In a "toy" MC measurement of the BLNY parameter g2, we study both the 
sensitivity to systematic uncertainties and the statistical sensitivity of qt, cil, 
Qt and 6^- Figure 5 shows the pythia (re-weighted to ResBos) prediction of 
the normalized distribution for two different values of g2- The ax distribu- 
tion is clearly sensitive to the value of g2. 

Event samples are generated using ResBos, for fifteen g2 values from 0.54 to 
0.82 (distributed around the world average, g2 = 0.68 toxn GeV^) [2]. Using the 
re-weighting procedure described earlier, pythia "MC samples" are produced 
corresponding to each of the g2 values. An independent sample of 200k pythia 
events is re-weighted to the central g2 value and represents the experimental 
"pseudo-data" sample. 

For each of the 15 MC templates, the pseudo-data vs MC is calculated. The 
function: y = a{x — bY + c is fitted to the as a function of g2, and a best fit 
g2 is determined as 6 ± a^^^"^ where the uncertainty is statistical (Ax^ = ±1). 

This "toy" measurement of g2 serves as an example analysis at an experiment. 
Free parameters in any NP model which affect the Qt could in principle be 
measured using this strategy. 
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7.1 Systematic Uncertainties 



We study the systematic uncertainties on a fit to g2. The following systematic 
variations are carried out. These are typical of the experimental uncertainties 
expected at a hadron colhder ahhough the size of the variations are chosen to 
give reasonable shifts in the fitted §2- 

• In the MC, the smearing constant is varied by a factor ±1% 
around a central value of 0.003 GeV~^. 

• As a test of sensitivity to mis-measurement of the event selection efficiency 
dependencies on, |?7|, px and isolation, events in MC are given ±10% more 
weight for each muon with; 1 < \ri\ < 2, 15 < p-r < 20 GeV, or 1 < fiso < 
2.5 GeV. 

• In the MC. the smearing constant is varied by ±10%. 

• As a test of sensitivity to mis-modeling of final state radiation (FSR), events 
in MC are given ±10% more weight if the difference between the generator 
level Z mass and the generator level di-lepton mass is greater than 1 GeV. 

Table 1 

The shifts (%) in the fitted g2 lor systematic variations in the AIC. 









Qt 


bx 


Pt smearing ±1% 


-0.02 


±11 


^2.8 


^0.03 


PT ±10% 


TO.2 


±3.8 


±1.4 


TO.2 


fiso ±10% 


TO.3 


±1.7 


TO.59 


TO.3 


\v\ ±10% 


±0.4 


±11 


±3.2 


±0.4 


(p smearing ±10% 


-0.01 


0.00 


0.00 


-0.02 


FSR ±10% 


TO.4 


TO.96 


TO.52 


TO.4 



Table 1 shows the shifts in the fitted g2 for the ±ve and -ve systematic vari- 
ations. For each of the variations, the shift in the fitted g2 is substantially 
smaller using qt or Bt than Qt, and larger using a^,. In fact the only variation 
to which Qt or bx are more sensitive is the smearing, but for any realistic 
detector resolution, the shift in g2 is negligible. 

7.2 Statistical Sensitivity 

Simply discarding information from is not optimal in terms of the statistical 
sensitivity to the shape of the Qt distribution (and thus also the value of 
g2). As well as the basic observables, Qt, aT and ql, the following ideas are 
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proposed as possible optimised combinations of ar and ai- 

• A weighted quadrature sum of Ot and Ol with more weight (w) given to ar'- 



• A 2D fit to -j^^. 
Table 2 

The binning and la statistical uncertainties for each of the g2 hts, for a range of 
5{1/pt). 





Qt ar aL QU^^ = 5) ^j£^ br 


nbins 
range (GeV) 


20 20 20 20 20(aT)x20(aL) 20 
0-30 0-20 0-20 0-30 0-20(aT), 0-20(0^) 0-20 


5{1/pt) (GeV-i) 


la statistical uncertainty (%) 


0.0000 
0.0003 
0.0010 
0.0030 


1.4 2.2 2.4 1.9 1.4 2.2 
1.4 2.2 2.5 1.9 1.4 2.2 
1.8 2.2 3.6 1.9 1.6 2.2 
3.1 2.3 8.9 2.2 2.1 2.2 



The binning and la statistical uncertainties for each of the fits is presented 
in Table 2. Since the relative statistical sensitivity of each of the observables 
depends on the detector resolution, the fit is performed for different widths of 
Gaussian smearing. All of the fits use equal width bins. For "perfect" resolu- 
tion, a fit to Qt and the 2D fit have comparable statistical uncertainties. A 
fit to Qt or ql alone is less statistically sensitive, as each contains information 
on only one component of the Qt- 

At larger S{1/pt), the sensitivity of ol is completely washed out by the smear- 
ing. For S{1/pt) — 0.003 GeV~^, ot alone gives better statistical precision than 
Qt- Over the resolution range covered, Qt{w = 5) gives slightly better statis- 
tical precision than Qt alone. Note that no attempt has been made to optimize 
the weight in Q^{w = 5) as a function of resolution, and w = 5 is a somewhat 
arbitrary choice. Clearly, for "perfect" resolution, w — 1 is the only sensible 
choice. Once experimental systematic uncertainties are taken into account, qt 
will need to be given more weight in any optimal combination of ot and ql- 



7.3 Sensitivity to Small-x Broadening 



The NP form factor (Eq. 1) required an alteration, to describe deep inelastic 
scattering (DIS) data involving initial state partons with x < 10~^ [13]. Berge 
et al. [14] suggested that if this (so-called small-x broadening) was observed 
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at the Tevatron in an exclusive high sample of Z bosons, a broader Higgs 
(and VT, Z) boson transverse momentum distribution could be expected at the 
LHC. The D0 Run II data on Qt [4] disfavoured this modification, although 
given the large uncertainties, the data was not particularly sensitive to such 
broadening. Even without taking into account the reduced systematic uncer- 
tainties, the optimised fits (Qt? (Per / daTdai) described earlier would be more 
sensitive to such effects. Once systematic uncertainties are taken into account, 
the sensitivity is further enhanced relative to the Qt- 

7.4 Azimuthal Correlation Between the Recoil and the Leptonic Decay 

The example g2 measurement presented would be sensitive to the descrip- 
tion by the MC event generator(s) of the azimuthal correlation between the 
hadronic recoil and the leptonic decay. A measurement of d^a / daxdai could 
be used to study this correlation. 



8 Conclusions 

Using MC simulations we demonstate the potential benefits of decomposing 
the Qt into two orthogonal components, ax and a^. A measurement of the 
ax distribution would be substantially less sensitive to two of the dominant 
experimental systematic uncertainties reported in previous Qt measurements: 
lepton Pt 01 Et mis-measurement, and the Qt dependence of the event selec- 
tion efficiencies. A slightly different variable 6t is demonstrated to be similarly 
insensitive to these uncertainties. An optimal combination of ot and ol, giving 
more weight to ot gives, for any realistic detector resolution, better statistical 
sensitivity to the shape of the region of low Qt- Two possibilities are proposed: 
a weighted quadrature sum of and ai or a fit to d^a/daTdaL. 

A measurement of d^a/daTda^ could potentially probe the azimuthal corre- 
lation between the Z boson decay axis and the hadronic recoil. The partial 
differential cross sections; d'^a /daTdQ, d^a/daxdy and d^a / daxdQdy could 
also be measured, taking advantage of the reduced systematic uncertainties 
on Qt compared with similar measurements of differential cross sections with 
respect to Qt- Such distributions would be sensitive to small-x broadening 
effects. 
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